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In this paper we study how to attack through different techniques a perfect fluid Bianchi I model 
with variable G,and A. These tactics are: Lie groups method (LM), imposing a particular symmetry, 
self-similarity (SS), matter collineations (MC). and kinematical self-similarity (KSS). We compare 
both tactics since they are quite similar (symmetry principles). We arrive to the conclusion that 
the LM is too restrictive and brings us to get only the flat FRW solution with G = const, and 
A = 0. The SS, MC and KSS approaches bring us to obtain the following solution: G is a decreasing 
time function and A ~ t~''^, with A < 0, while the exponents of the scale factor must satisfy the 
conditions X]i=i = 1 and X]i=i Q^? < li ^ 1) i relaxing in this way the Kasner conditions. 

PACS numbers: 98.80.Hw, 04.20.Jb, 02.20.Hj, 06.20.Jr 



Contents 



I. Introduction. 2 

II. The Model. 3 

A. The crucial equations. 4 

III. Curvature behavior. 6 

IV. Lie Method. 7 

A. Scale symmetry. 8 

B. Exponential behavior. 9 

C. Solution with the full symmetry. 10 

V. Self-similar solution. 11 

VI. Matter collineations. 15 

A. The usual matter collineation equations. 15 

B. Modified matter collineations equations 16 

C. The complete modification of MC equations. 17 

VII. Kinematic Self-similarity. 17 

VIII. Conclusions. 21 

Acknowledgments 22 

References 22 

A. Lie Method. 23 

1. Scaling symmetry 24 

B. On self-similar solutions. 26 

1. "Constants" constant. 26 

2. G-variable. 27 

3. A-variable. 28 

4. G&A-variable. 29 



'Electronic address: abelcal@ciccp.es 



2 



I. INTRODUCTION. 

Ever since Dirac first considered the possibility of a G variable (see [1]), there have been numerous modifications 
of general relativity to allow for a variable G, nevertheless these theories have not gained wide acceptance. However, 
recently (see [2]-[ll]) a modification has been proposed treating G and A as non-constants coupling scalars. So it is 
considered G and A as coupling scalars within the Einstein equations, Rij — ^gij = GTij — Ag^j, while the other 
symbols have their usual meaning and hence the principle of equivalence then demands that only gij and not G and 
A must enter the equation of motion of particles and photons, in this way the usual conservation law, divT = 0, 
holds. Taking the divergence of the Einstein equations and using the Bianchi identities we obtain the an equation that 
controls the variation of G and A. These are the modified field equations that allow to take into account a variable G 
and A. Nevertheless this approach has some drawbacks, for example, it cannot derived from a Hamiltonian, although 
there are several advantages in the approach. 

There are many publications devoted to study the variation of G and A in the framework of flat FRW symmetries 
(see for example [2]-[ri]) and all this works have been extended to more complicated geometries, like for example 
Bianchi I models, which represent the simplest generalization of the flat FRW models (see for example [12]-[15]. in 
the context of perfect fluids and [10]-[18] in the context of viscous fluids). Bianchi I models are important in the study 
of anisotropics. 

But in our opinion, the problem arises when one try to solve the resulting field equations (FE). It seems that it is 
unavoidable to make simplifying hypotheses, or to impose ad hoc some particular behaviour for some of the quantities 
of the model, in order to obtain a exact solution to the FE Such simplifying hypothesis are made for mathematical 
reason (in order to reduce the number of unknowns) although are justified form the physical point of view. Usually 
such assumptions or simplifying hypothesis follow a power law, for example, the quantity X follows a power law i.e. 
X = Xot", where Xq is an appropriate dimensional constant, t is the cosmic time (for example) and a G E (usually 
a £ Q, but this other question), and depending on the nature of the quantity X, a will be positive or negative. 
Actually we think that although all these simplifying hypothesis are correct or at least bring us to obtain correct 
results, it is not necessary to do that, since they may be deduced from symmetry principles in such a way that one may 
justify (deduce) them from a correct mathematical principle, and usually all these approaches have physical meaning. 

Therefore the main goal of this paper is to develop and compare some well known tactics (approaches) in order to 
study and find exact solutions for a perfect fluid Bianchi I models with variable G and A, but trying to make the 
lowest number of assumptions or neither. We will try to show that with these approaches all the usual simplifying 
hypotheses may be deduced from a correct mathematical principle and how the useful are each tactic, i.e. to show 
the advantages and disadvantages of each approach. We have started studying this class of models because, as we 
have mentioned above, there are many well known exact solutions so we will be able to compare the useful of our 
approach. In a forthcoming paper we extend this study to a very complicated model (from the mathematical point 
of view) which is a perfect fluid Bianchi I within the framework of a variable speed of light (VSL) where we have to 
study, for example a third order ODE with four unknowns, and we need to be sure that theses tactics, the exposed 
ones here, work well. Hence in this paper we are going to study a Bianchi I model with variable G and A through the 
Lie group method (LM), studying the symmetries of the resulting ODE, and through the self-similarity (SS), matter 
coUineations (MC) as well as kinematical self-similarity (KSS) hypothesis. 

The paper is divided in the following sections: In section two we outline the main ingredients of the model as well 
as the FE (under the condition divT = 0). In order to apply the (LM) we need to deduce a ODE. For this purpose we 
have followed model proposed by Kalligas, Wesson and Everitt, [13], but taking into account some little differences. 
In section three are calculated all the curvature tensors, Weyl etc... as well as their invariants, i.e Kretschmann scalars 
etc.... In section four it is studied through the Lie group tactic a third order differential equation with two unknowns. 
We seek the possible forms that may take G in other to make integrable the ODE. In this way we flnd that there are 
three possibilities, but the question here, is that all the studied solutions are unphysical (in the sense that the shear 
vanish i.e. ct = 0) or trivial. Actually this is a very surprising result, since the third solution is very similar to the 
obtained one in reference [13], but if one follows all the calculations until the end, then arrive to the conclusion that 
such solution is the flat FRW one with G constant and with a vanishing A. Since in the case of a c— var the third 
order ODE is quite complicate, in appendix A, we study the resulting second order ODE, which is simpler than the 
third one but has the drawback of having three unknowns. In this case is really difficult to arrive to some conclusions 
since the obtained solutions depend of many integrating constants although the solutions have or follow the same 
behavior, in order of magnitude, as in the previous case. 

In section five, we study the model under the self-similarity hypothesis, in this case, the obtained solution are 
similar (in order of magnitude) to the obtained one in the above section (LM with scaling symmetries). Nevertheless 
in this case we are able to obtain a solution with shear non- vanishing i.e. ct 7^ 0, where G is a decreasing time function 
and A « with A < 0, while the exponents of the scale factor must satisfy the conditions cti = 1 and 

X]i=i Q^f Ij '^^ ^ 1) I relaxing in this way the Kasner conditions.. Since the model is SS, then, in section six. 
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we study the model studying the matter coUineations (MC). In this occasion we need to reformulate or reinterpret 
the MC equations in order to get information on the behavior of G and A, arriving to the same conclusions as in 
the above section. In the last section we reproduce the same tactic but this time under the KSS hypothesis, in this 
occasion we get a non-singular solution and with the same behavior for the main quantities as the obtained one in 
the above sections. We end this section discussing the Kasner like solutions. We end with a brief conclusions. 

We have add another appendix, appendix B, where we review a set of solutions for Bianchi I models with time 
varying constants but under the hypothesis divT ^ 0, as well as, the standard one, which is the model where G = const. 
and A = 0. We arrive to similar conclusions as the obtained ones in the paper. 



II. THE MODEL. 



Throughout the paper M will denote the usual smooth (connected, Hausdorff, 4-dimensional) spacetime manifold 
with smooth Lorentz metric g of signature (—,+,+,+). Thus M is paracompact. A comma, semi-colon and the 
symbol C denote the usual partial, covariant and Lie derivative, respectively, the covariant derivative being with 
respect to the Lcvi-Civita connection on M derived from g. The associated Ricci and stress-energy tensors will 
be denoted in component form by Rij{^ R'^jcd) and Tij respectively. A diagonal Bianchi I space-time is a spatially 
homogeneous space-time which admits an abelian group of isometrics G3, acting on spacelike hypersurfaces, generated 
by the spacelike KVs ^1 =8x^2 = dy,£,3 = dz- In synchronous co-ordinates the metric is: 



(1) 



where the metric functions Ai{t), A2{t), ^3(i) are functions of the time co-ordinate only (Greek indices take the space 
values 1, 2, 3 and Latin indices the space-time values 0, 1, 2, 3). In this paper we are interested only in proper diagonal 
Bianchi I space-times (which in the following will be referred for convenience simply as Bianchi I space-times), hence 
all metric functions are assumed to be different and the dimension of the group of isometrics acting on the spacelike 
hypersurfaces is three. Therefore we consider the Bianchi type I metric as 



ds-" 



c^df 



+ X^{t)dx^ + Y^{t)dy^ + Z^{t)dz 



(2) 



see for example ([19] -[2 5]). 

For a perfect fluid with energy-momentum tensor: 



Tij = (p + p) UiUj + pgij, 



(3) 



where we are assuming an equation of state p = up, [lo = const.). Note that here we have preferred to assume this 
equation of state but as we will show in the following sections this equation may be deduced from the symmetries 
principles as for example the self-similar one. The 4— velocity is defined as follows 



u = ( i, 0,0,0 



UjU^ = — 1. 



(4) 



The time derivatives of G and A are related by the Bianchi identities 



Ri 



:R9i. 



( 



SttG 



V c4 



Tij - hgi 



(5) 



in this case this equation reads: 



SttG 



p + p(l + ^) 




-A 



8^ . 

-^Gp, 



(6) 



in our case we obtain (taking into account the additional condition (divT = 0)) : 



, X Y Z\ 



A = — jGp, 



(7) 
(8) 
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in the appendix B we will study how to obtain different cosmological models, under the self-similarity hypothesis, 
with one or several time varying "constants" but without the assumed condition div T = 0. 
Therefore the resulting field equations are: 

XY XZ ZY SirG . . 

Y Z ZY 8ttG ^ 2 

Y + Z + ZY=-—^' + '''^ (10) 



X Z XZ 8ttG ^ ^ 

x + z + xz--— + (11) 

X Y XY _ SttG 
X^Y^ XY ^ ^ 



ujp + Ac, (12) 



p(l+^)(^ + I + |) =0, (13) 



X Y Z 



A=-^Gp. (14) 



Now, if we define 



X Y Z\ R o d { l\ 

i/=(^-^ + ^ + ^j^3- and R^^XYZ^ '^^^Uj-l' ^'^^ 

it is regained eq.(13) as the usual conservation equation i.e. 

p + p{l+uj)H = Q. (16) 
Since we have defined the 4-velocity by eq. (4) then the expansion 6 is defined as follows: 





Y 




1 




'y' 


-D 


c 



^ <^ (^--A- + ^ + -\=zH, (17) 

and therefore the acceleration is: 

at = Ui-jU^ , (18) 
in this case a = 0, while the shear is defined as follows: 

CTij = ^ {Ui-ah'^ + Uj-aK) - \^9hij , (19) 

or equivalently 





X 










z 


X 










z 



-' = 7^11^-^) + (^-^1 + I- (21) 



A. The crucial equations. 

In this section we would like to obtain a ODE which allows us to study all the field equations through the Lie 
method. For this purpose we are following closely the paper by Kalligas et al (see [13]). From eqs. (10-12) and taking 
into account cq. (9), we get: 

f + 7 + | = -5(l + 3'^)Gp + Ac^ (22) 
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Now, taking into account cq. (13), squaring it we get 



-) =(l + c^fi^^ (23) 



and taking into account (9) we get 



The time derivative ^ from eq. (13) can now be expressed in terms of G, A and p only by using eqs. (22) and (25), 
it is obtained: 

^ 94-,, I n ^9■w 

-Gp{l-uj) + 3Ac^. (26) 



2 / . \ 2 / . \ 2 



pJ (l+u;) 



Now if we rewrite this equation in an appropriate way we get the foUowing one that we wiU also study through the 
Lie method (see appendix A)., i.e. 

P = A^i — + K2Gp'- - KsAp, (27) 
P 

where 

Ki = jT^, K2 = ^ i^3-3 1+c.c2. 28 

+ U)) 

as it is observed if cj = 1, then K2 = 0, so eq.(27) collapses to a very simple ODE with only two unknowns. 

In order to reduce the number of unknowns, now, on differentiating eq. (26) and taking into account eq. (14) we 
arrive to the following equation 



. 3 

'p = Kip^ -K2^+ KiGp^ - KiGpp, (29) 



where 



5 + 3^ 4 + 2^ 12^(1 + ^)^ 12^(1-0^^) 

^^ = TTT7' ^'^ = TTT7' = 7^ ' ' ^^^^ 

we are supposing that uj ^ —\. This equation i.e. eq. (29) will be studied exhaustively through the Lie method in 
section IV. 

As it is observed if = 1 (ultra-stiff matter) eq. (29) collapses to 

■P = Apt -3^ + '^Gp\ (31) 

The shear is calculated as follows. If we take into account the definition of the shear i.e. eq. (21) as well as eqs. 
(9), (13) and (24) then we arrive to the following expression 



- 1^'"' - (« V + ^} - (y - - " 

note that lu ^ —1. We are only interested in solutions which verify a ^ 0. 

Remark 1 In reference ([1-^]) the authors define the shear in a different way (see eq. I4 of reference ([1 '']))■ They 
define the anisotropy energy as 

S^Ga = Q(^[j^ - (^Stt^P + A^^ , (33) 
which is similar to our definition of shear (compare with eq. 32). 



III. CURVATURE BEHAVIOR. 

The full contraction of the Riemann tensor (see for example [26]-[29]) 

Ii RijkiR^-'^^ , 



h 




The full contraction of the Ricci tensor 




and the scalar curvature R is: 



X'Y' Y'Z' X'Z' 



XY 



YZ 



XZ 



We calculate the Weyl tensor as well as the magnetic and the electric part of the Weyl tensor. 
The non-zero components of the Weyl tensor are: 

X" X' (Y' 



C1212 
C1313 

Cl414 
C2323 
(^2424 
^3434 



6 
6 

ZP_ 
6 

6^ 

6^ 

y2^2 



X 
Y 
Z 



X 

Y 



2— + — — + — + — 



6c2 



Y 

'x 

'x 
z^ 

XL 
Y 

X 



z^ 

z 
z^ 

" Y 

Y 

)C_ 

Y 

}C_ 

Y 
r' 

Y 



Y 
x^ 
'Y 
x^ 
'Y 



-2^ 
z 

Y" 



Z" 


-4 






^Y 




f) 




Z" 


-2^ 






^Y 


X 


5) 




Y" 










X 


y) 




X" 


Y" 






-Y^ 


Y' ^ 


X 




X" 


Z" 


X' 




-Y^ 


Y ^ 


^Y 


1) 


Y" 


Z" 


2XL 






Y ^ 


Y 


§) 



where X' := X. 

The non-zero components of the electric part of the Weyl tensor are: 



^^22 — 



33 



Eaa — 



X^ 

6c2 

y2 

6c^ 
Z^ 



-2- 



-2- 



r" 
) 



X^ fY^ 

Y [y 



XL 
Y 

XL 
Y 



Y 

XL 
Y 



XL 
Y 

XL 
Y 

XL 



+ 



XL 
Y' 

XL 
YX 

XL. 
'Y 



Z" Y' Z' 

Y ^ '^YY 



XL 
Y 

Y" 



iXLXL 
'YY 

X' Y' 



Y '^'X Y 
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The magnetic part of the Weyl tensor vanish 



The Weyl scalar is defined as: 



J3 := C'^'d^a. 



hcdi 



(48) 
(49) 



4 

3^ 



as it is observed is also defined as follows: 




Z" 


Y" 


~z " 


' T 








~z) 



(50) 



h = h- 2/2 + 

this definition is only valid when n = 4. 

The gravitational entropy is defined as follows (see [27]- [28]): 



P2 



h 
I2 



h - 2/2 



1 d2 
3-^ 



h 
h 



(51) 



(52) 



In this section wc will study eq. (29) i.e. 



IV. LIE METHOD. 



(53) 



where the (A'i)^^]^ are given by eq. (30), through the Lie group method, i.e. in particular we seek the forms of G{t) 
for which our field equations admit symmetries i.e. are integrable (see for example [3()]-[34]). 
Following the standard procedure we need to solve the next system of PDEs 





= 0, 


(54) 




= 0, 


(55) 


Kipri - Kip^rjp ~ 9p^^tp + ^P^Vpp 


= 0, 


(56) 


-KlP^T]t + Sp^TJip - ip^^tt 


= 0, 


(57) 


Kip'^^pp + K2p^p - p^(.ppp 


= 0, 


(58) 


2K2PVP - Kip^r]pp + 2Kip^^tp - " ip^^tpp + P^Vppp 


= 0, 


(59) 


3K2PVt - 2Kip'7jtp + SK^p^G^p + Kip^^tt - ^pHttp + ip^ltpp 


-0, 


(60) 


Kip^G'C + Kip^Gr) + 2Kip''G^t - Kip^rjtt + ip^'vup - P^^ttt - ^K^P^'G'^p 


= 0, 


(61) 


Kip^Grtt + p\tt + K^P^G'ijp - SK^p^G'^t ~ K^p^G''^ - 2K3p^G'v 


= 0, 


(62) 



To study all the possible solutions for this system of PDE may be very tedious, for this reason we impose the 
symmetry X ~ {at + e) dt + bpdp, i.e. = at + e, rj = bp, in such a way that we get the following restrictions for G{t). 
From eq. (61) we get 



G^ 
G 



6+ 2a 

at + e' 



(63) 
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while from eq. (62) it is obtained: 

G" _ 3a + 6 
~ ~ at + e ' 



(64) 



where a,6, e S M. Note that [a] = [6] = 1, i.e. they are dimcnsionlcss constants but [e] = T, with respect to a 
dimensional base B = {L, M, T} . So we are going to study the different solutions that we may obtain in function of 
the constants a, 6, e. 



A. Scale symmetry. 

Making e = i.e. considering only (^ = at, r/ = 6p) , we have to integrate eqs. (63 and 64), so 

=^ G = Go(t)-('+»), (65) 
=^ G = G2 + Go(i)"('+"), (66) 



G' _ b + 2a 
"G ^ 

G" 3a + & 



G' at 
therefore we get 



G = Go (0"('+") , (67) 



where we assume that Go > 0. 

The invariant solution for the energy density is: 

bdt dp ,,/„ 

— = — ^ P = Po^'/", (68) 
at p 

and for physical reasons we impose the condition, ab < then & < 0. We have considered only the invariant solution 
since as we already know, the general one usually lacks of any physical meaning as in this case (see for example [-34]). 
Furthermore, as we will see in section V, this spacetime is self-similar, this means that all the quantities must follow 
a power law as in this case (see for example [38]-[39]). 

If we make that this solution verifies eq. (53) with G{t) given by eq. (67), we find the value of constant po, so 

c^b 

Pa = n, (69) 

Gol27ra(l + w)' 

with the only restriction uj —1. Note that ab < 0, so po > 0. 

Remark 2 As we can see, it is verified the relationship Gp = t^^, i.e. the Mach relationship for the inertia. 

Therefore, at this time we have the following behavior for G{t) 

decreasing if b/a G (—2, 0) , 
G{t) = Got-^^'^-\ G w { constant if b/a = -2, . (70) 

growing if b/a < —2 

With these solutions we find the behavior of the rest of the quantities. From eq. 

TT ■ Stt /„ 5\ ^ _n . . _9 . 1 b ^ b 



I.e. 



A . --^Gp = ^ ^2 + - j poGot-^ =^ A . Aot-, Ao = ^J^^- [2 + - ) , (71) 



negative if 5/a e (—2, 0) , 
A(t)==Aot"^, Ao « <( vanish if 6/a -2, , (72) 

positive if b/a < —2 

therefore, if Ao > 0, then G is a growing time function but if Aq < 0, then G is decreasing. 
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H = -|.7T^I9- . = = (73) 



Whit regard to H we find that from eq. (15) 

a{l+uj)) i ' " ^ dt \H 
so 

R = i?op"^/^^^+'^^ = i?oi-''/^"(^+'^) , Z = i?ot~''/"^^+"^ . (74) 

If we assume that the functions (X, F, Z) foUow a power law (i.e. X ~ Xot"^) then we get the following result 

Kt- = R,t-^l<'^-\ =^ (75) 

i 

where we may assume that [ai) > 0,Vi and (a^ 7^ aj) although (a^ aj) when t — > 00, with i ^ j, but we have not 
more information about this behavior. 
The shear is calculated as follows. 



2 



3c2(l+cj)' (p) V^c4^^^) 3c^{l+LjfW^'^a aT^aj^^ ^^^^ 



therefore the shear vanish. In this way. this is an unphysical solution since we are only interested in solutions which 
verify the condition a 0. Actually, we have obtained the flat FRW solution as we will show at the end of this section. 

Nevertheless if we follow the hypothesis of power law for the scale factors {X, Y, Z) , then we get from the definition 
of the shear (see eq. (20)) that it follows 

note that a — > when (ai —faj). Note that at least one expect this behavior but as we will see in the bellow 
subsection (full symmetry) actually a = 0, since ai = aj, (the fiat FRW solution) so this is an unexpected result. 

B. Exponential behavior. 

Making a = 0, we have {£, = e, r] = bp) , so integrating eqs. (63-64) it yields 

% = --,=^G = GocM--t), (78) 
G e e 



in this way we obtain 



%^--,=^G = C2 + G„cxp{~-t), (79) 
G e e 



G = Goexp(--t), (80) 
e 



where we assume that Go > 

The invariant solution for the energy density is: 

bdt dp ,h . 

— = — =^ p = poexp(-t) (81) 

e p e 

with the restriction, eb < with 6 < 0, from physical considerations. In order to calculate the value of constant 
po, this solution must verifies eq. (53) finding in this way that constant po, vanish i.e. po = 0. Note that with the 
obtained behavior for G and p, such functions must verify the field eq. (14). 
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C. Solution with the full symmetry. 

In this case we have the full symmetry i.e. = at + e, rj = bp) , so, by integrating the constrains we find that 

^' .G = Go(a^ + e)-(2-+'')/^ (82) 



G at + e 

G" _ 3a + 6 

^ ^ at + e 



G = Go(a^ + e)-(2'^+'')/^ (83) 



with [e] = T, that is, we may interpret it as a characteristic time, hke in the ad hoc sohition given by KaUigas et al 
(see [13]). 

The invariant solution for the energy density is 

p==po(at + £)"/", (84) 

where we need to impose the physical constrain such that a& < then & < 0, in such a way that it must verifies eq. 
(53) with G{t) given by eq. (83), finding the value of the numerical constant poi so 

c% , , 

Po = 7, (85) 

12naGo{l + Lu) 

we assume that uj ^ —1. As it is observed this is a nonsingular solution since when t ^ if e 7^ 0, then p ^ 00. 
Therefore we have the following behavior for G(t) : 

growing if 5/a G (—2, 0) , 
G{t) = Go(at + e)-(2a+fc)/a^ q _ ^ constant if b/a = -2, . (86) 

decreasing if b/a < —2 

Remark 3 Note once again that it is verified the following relationship: Gp » (at + e)^^. 
We go next to calculate the rest of the quantities. From (14) we get 

i.e. 

{negative if \b\ < 2a, 
vanish if |6| = 2a, . (88) 
positive if \b\ > 2a 

With regard to the quantity H, we find form eq. (15) that 

R = Rap-^/'^^^+^^ = Ro{at + e)-''/3a(i+^)^ (gg) 

and 

/ b \ I d f 1\ a(l+u;) 

hence 

XYZ ^ Ro{at + ey''/''^^+'^\ (91) 

so (following the same argument as above) the functions {X, Y, Z) follow a power law (i.e. X = Xo(at + e)""^, etc..) 
it is found that 

3 , 

K{at + er^Ro{at + e)-'/^^'+-\ =^ Y^a,=a = -^^^, (92) 
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where we may "assume" that (a^) > 0,Vi and (a^ ^ aj) although (ai — > aj) when t oo, and i ^ j. 
The shear has the foUowing behavior. 

^ -(8n^P + A)= Jb' (l - ^) +-) (at + e)-'^ ^0, (93) 

3c^^+uJf\pJ V J 3c2(l+cc>)H V aj' 

finding that 

6 (6 (2a2 - 1) + 2a) > b {2a^ ^ l) + 2a < 0, 4=> \b {20^ - l)\ > 2a =^ ( " ^ ' (^"^^ 

since CTq must be positive. 

At this point it seems that we have found a physical solution that depends on the value of constants a and b. But, 
how to calculate the value of constants (at)^^^?. 

Now we will try to calculate the possible values for the constants (ai)'^^^ , simply all these results must satisfy the 
field eqs., so they are solution of the following system: 

1/^2 2a \ 

aia2 + aias + a2a3 = - [ a + -— — r , (95) 
I f 2 2a{2Lo + I] 



a2 {a2 - 1) + as (as - 1) + q;3Q;2 = tt " 77-, ^ — , (96) 

\ [1 + ijj) 

1(2 2a (2^ + 1) 



ai (ai - 1) + as (03 - 1) + a^ai = 77 a 7— — r — , (97) 

I f ■) 2a{2uj + l)\ 

ct2 ("2 - 1) + ai [ai - 1) + aia2 = 77 a , (98) 



6V (1 + ^^) 

and this system has only two solutions, the trivial one {ai = 0:2 = ^3 = 0) and the flat FRW one i.e. 

2 

Qfi = Qfa = ^3 = ^7^— — (99) 
3 (1 + w) 

which is incredible, this means, that from eq. (92) we obtain b ~ —2 and a = 1, but with these values for constants 
(a, 6) , G is a true "constant" , A vanishes and (7 = 0. This solution was obtained by Einstcin&de Sitter ([()■!]) in 1932 
for u — 0, and later by Harrison ([64]) 

We think that the followed method is too restrictive and for this reason we obtain only this solution. 

Therefore this method brings us to obtain only the classical flat FRW solution with G constant and A equal nought. 



V. SELF-SIMILAR SOLUTION. 



In general relativity, the term self-similarity can be used in two ways. One is for the properties of spacetimes, the 
other is for the properties of matter fields. These are not equivalent in general. The self-similarity in general relativity 
was defined for the first time by Cahill and Taub (see [35], and for general reviews [36]-[42]). Self-similarity is defined 
by the existence of a homothetic vector V in the spacetime, which satisfies 

Cvgij = 2agij, (100) 

where Qij is the metric tensor, Cy denotes Lie differentiation along V and a is a constant. This is a special type 
of conformal Killing vectors. This self-similarity is called homothcty. If a ^ 0, then it can be set to be unity by a 
constant rescaling of V. If a = 0, i.e. Cygij = 0, then y is a Killing vector. 

Homothety is a purely geometric property of spacetime so that the physical quantity does not necessarily exhibit 
self-similarity such as CyZ = dZ, where d is a constant and Z is, for example, the pressure, the energy density and 
so on. From equation (100) it follows that 

CvR'jki^O, (101) 

and hence 



CyRij = 0, CyGij = 0. 



(102) 
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A vector field V that satisfies the above equations is called a curvature coUineation, a Ricci collineation and a matter 
collineation, respectively. It is noted that such equations do not necessarily mean that y is a homothctic vector. We 
consider the Einstein equations 

Gij^SnGTij, (103) 

where T!y is the energy-momentum tensor. 

If the spacetime is homothetic, the energy-momentum tensor of the matter fields must satisfy 

CvT,, = 0, (104) 

through equations (103) and (102). For a perfect fluid case, the energy-momentum tensor takes the form of eq. (3) 
i.e. Tij = {p + p)uiUj +pgij , where p and p are the pressure and the energy density, respectively. Then, equations (100) 
and (104) result in 

Cyu' = —au\ CvP = --2ap, Lyp = — 2ap. (105) 

As shown above, for a perfect fluid, the self-similarity of the spacetime and that of the physical quantity coincide. 
However, this fact does not necessarily hold for more general matter fields. Thus the self-similar variables can be 
determined from dimensional considerations in the case of homothety. Therefore, we can conclude homothety as the 
general relativistic analogue of complete similarity. 

From the constraints (105), we can show that if we consider the barotropic equation of state, i.e., p — f{p), then 
the equation of state must have the form p = cup, where is a constant. This class of equations of state contains 
a stiff fluid {lo = 1) as special cases, whiting this theoretical framework. There are many papers devoted to study 
Bianchi I models (in different context) assuming the hypothesis of self-similarity (see for example [43]-[44]) but here, 
we would like to try to show how taking into account this class of hypothesis one is able to find exact solutions to the 
field equations within the framework of the time varying constants. 

The homothetic equations are given by eq. (100) so it is a straightforward task to find the homothctic vector field, 
where in this case is as follows: 

/ v\ I \r\ I 'vx 

(106) 



XX + tXX-t[X] )a; = 0, (107) 




iff the following ODE is satisfied 



and so on with respect to (F, y) and (Z, z) . 

As it is observed from eq. (107) if we simplify this ODE it is obtained the following one: 

^ = — tFi = const. X = Xot"^ , (108) 
Hi t 

with Xq and a G M, etc....with respect to the others scale factors. Note that ' '■— ■= dot. i.e. X' = X. 
Therefore, we have obtained the following behavior for the scale factors: 

X = Xot"^ , Y = Fot"' , Z = Zot"-^ , (109) 

with Xq, Yq, Zq are integrating constants and {cti)^^i G K. In this way we find that 



^= I >J"» I 7 = ("i+"2 + a3)y = (110) 



and hence 



p = (Ill) 
From the field equations (9) and (14) we get that 

(^(j,^-{u+i)a _ _^ l)aGr , (112) 



^, _ A 2 8npo 



c2 i3 
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where A = X^iT^ i aiCtj = aiQ;2 + ctac^i + oi2Ci-i, a = ai = ai + a2 + 0:3, therefore 



■2 A 



G=-— 7 ^i(-+i)"-2 (^^3) 



as we can see it is verified the relationship Gp « t ^ as it is expected. 
The behavior of G is the following one: 



growing if + 1)q! > 2 
G = G'ot'"+^^""^ =^ { constant if + l)a = 2 , (114) 

decreasing if (w + 1)q! < 2 



with Go > 0. 

Now, we go next to calculate the quantity A, from eq. (112) we get 



A = Aot-2, Ao = 4 ( 1 - r^T^ ) ' (115) 

' [oj + l)a ' 



in this way it is observed that 



Ao > <^ {lo + 1)q; > 2 
A = Ao^"^ <( Ao = ^ (cj + l)a = 2 . (116) 
Ao < {uj + l)a < 2 

The shear behaves (see eq. (20)) as follows: 

^' = ^ I I ^^0' (117) 



by hypothesis, since at this point we have not any information about the value of the constants {ai)^^i ■ 

In order to find the value of constants (0;^) , we make that them verify the field eqs. so in this case we get the 
following system of eqs.: 



a2 (a2 - 1) + as (as - 1) + a3a2 = A y — - — j , (US) 

ai (ai 1) + as (as - 1) + aaai = A ^ ^ , (119) 

a2 (a2 - 1) + ai (ai - 1) + aia2 = A [- — -) , (120) 



where A = Q;ia2 + a^ai + a2as, and a = ai + a2 + as. 

So we have the following solutions for this system of equations: 

"1 = "2 = as, (121) 
ai = 1 - a2 - as, (122) 

as it is observed solution (121) is not interesting for us, since it is again the flat FRW solution. Only the second 
solution has physical meaning (in this framework). Nevertheless we have found that this solution only verifies the 
first of the condition of the Kasner like solutions i.e. 

a = ^a, = 1, (123) 

while the second condition 

it is not verified (see [Gl] and [G2]). In this case we find that it is verified the condition 

E"?<1- (125) 
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Therefore we have found a similar behavior as the obtained one in ([43], we say similar because there the authors 
only study standard models i.e. models where the "constants" are true constants, in fact A = 0), except than here 
this result is valid for all equation of state i.e. Voi G ( — 1, 1) . Nevertheless in reference ([44]), the authors claim that 
the solution must verify both conditions, i.e. ^o^i = 1 = ^af. See the end of the kinematical self-similar solution 
for a comment on this class of solutions as well as the appendix (B). 

Therefore we have obtained the following behavior for the main quantities: 

H=j, =^ q = 0, (126) 

so it is quite difficult to reconcile this model with the observational data. With regard to the energy density we find 
that 

p = Poi-^"+'\ (127) 

so, if w < — 1 =^ p is growing (phantom cosmologies), for the rest of the values of uj, i.e. w G (— 1, 1], p is a decreasing 
function on time. 

G has now a more restrictive behavior since (we are supposing that lu G ( — 1, 1], i.e. w > 1, is forbidden) 

„ _ „ j^+i)-2 , J constant if (w = 1) , . 

^~^ot , =^ \ decreasing Vc^G (-1,1) ' ^^^^> 

therefore we find that G is a decreasing function on time. The cosmological "constant" behaves as follows 

^-^"^-^ ^"-{t;S,""v:G(-i,i) ' (129) 

so we have found that A is a negative decreasing function on time. 

As we can see, this solution is quite similar to the obtained one through the Lie method with the scaling symmetry, 
at least in order of magnitude. 

With regard to the curvature behavior, we may see that 

^1 = ^2 = , (130) 

where K = K{ai) = const ^ 0, i.e. 

K = [3 {a\ + o?^ + 2a2a3 + 9a^a^ + 3 (a^ + aj^) - 6 (af + a^) + a^a^. (6 {a\ + ag) - 8 {a-i + as))] . (131) 
The non-zero components of the Weyl tensor are: 

C1212 = A'it-2("2+a3)^ = /^2t-'^l-"^\ Ci414 = K^r''^^- '^-\ 

^2323 =^4^"'"% C2424 - A-5^-2"^ C3434 = (132) 

where {Kif^^Y = K{ai) = const ^ 0, and taking into account a very famous result by Hall et al (see [47], and the 
next section VI) we may check that 

CvC%i = 0, (133) 

as it is shown in ([47]) if a vector field V G X(M), verifies the conditions CyCju = 0: ^^'^ CvTij ~ 0, then Cyg = "^g 
i.e. it is a homothetic vector field, but in this case we have arrived to the conclusion that Cyg = 2g 4=4' CyTij = 
(as we will see in the next section) and that it is also verified the relationship CyC'^ju = 0. 
The non-zero components of the electric part of the Weyl tensor are: 

E22 = A'it~2(a2+a3)^ ^ /^2<"^^^"""\ S44 = /l3<"^^^""=^\ (134) 

and the last invariant has the following behavior 

h = (135) 



with K given by 



k = -^\a\ + a\~ a2az + ia\oi\ + a\ + a\-2 {a\ + a\) + 0203 (2 {a\ -f 0:3) - {a2 + ^3))] , (136) 



3 

while the gravitational entropy behaves as 
note that P'^ = /3//2. 



p2 = const. ^ 0, (137) 
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VI. MATTER COLLINEATIONS. 



In recent years, much interest has been shown m the study of matter cohmeation (MCs) (see for example [45]-[53]. 
A vector field along which the Lie derivative of the energy-momentum tensor vanishes is called an MC, i.e. 

CvT.j = 0, (138) 

where V"^ is the symmetry or collineation vector. The MC equations, in component form, can be written as 

fcF^- + T^kV!] + Tk, = 0, (139) 

where the indices k run from to 3. Also, assuming the Einstein field equations, a vector generates an MC if 
CyGij =0. It is obvious that the symmetries of the metric tensor (isometrics) are also symmetries of the Einstein 
tensor Gij, but this is not necessarily the case for the symmetries of the Ricci tensor (Ricci collineations) which are 
not, in general, symmetries of the Einstein tensor. If F is a Killing vector (KV) (or a homothetic vector), then 
CyTij = 0, thus every isometry is also an MC but the converse is not true, in general. Notice that collineations can 
be proper (non-trivial) or improper (trivial). Proper MC is defined to be an MC which is not a KV, or a homothetic 
vector. 

Carol et al (see [4()]) and Hall et al.(see [47]) have noticed some important general results about the Lie algebra of 
MCs. 

Let M be a spacetime manifold. Then, generically, any vector field V on M which simultaneously satisfies CyTab = 
(<^ CyGab ~ 0) and CyC^^^ = is a homothetic vector field. 

If is a Killing vector (KV) (or a homothetic vector), then CyTab = 0, thus every isometry is also an MC but the 
converse is not true, in general. Notice that collineations can be proper (non-trivial) or improper (trivial). Proper 
MC is defined to be an MC which is not a KV, or a homothetic vector. 

Since the ST is SS then we already know that the SS vector field is also matter collineation i.e. we would like 
to explore how such symmetries allow us to obtain relationships between the quantities in such a way that it is not 
necessary to make any hypothesis to a solution to the field equations. In order to do that we need to modify the usual 
MC equations since with the usual one we are not able to obtain information about the behavior of G and A. 



A. The usual matter collineation equations. 

For a vector field V = {Vi {t, x, y, ^))^^i , the matter collineations reads: 

CyT.j = 0, (140) 

so outlining the equations and integrating them (it is a straightforward task) we obtain the following interesting 
relationships 

- = -2^, ^p^KoV,-^. (141) 



and 



in this way we find that 



^3 = -vi ( ^1 - ^) y> ^4 = -Vi (^1 - ^ ) ^> (143) 



therefore the obtained matter collineation (MC) vector field is: 



V = Vidt + VA^- Hi) xd,. +Vi(^- H2) ydy + Vi(^~H3] zd,, (144) 



Vi J \Vi 7 " v^i 

where as it is observed if Vi = t, then it is obtained the homothetic vector field (see eq. (106)) as we already know 
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For this reason wc may also check that the homothetic vector field verify the relationship as well 

LHoT^, = 0, (145) 

iff the following ODEs are satisfied: 

c2 {p't + 2p) = 0, {p't + 2p) = 0, (146) 

-px (^XX' + tXX" - t {X' f^ = 0, (147) 

as we can see the ODEs (146) {p't + 2p) = 0, and p't + 2p = brings us to obtain p = ujp » t~'^, as it is expected in 
a model without time varying constant and from the eq. XX' + tX" — t {X') i.e. Hi + tH'i = 0, that tHi = const. 
and so on with respect to the scale factors Y and Z i.e. tHi = const, i = 2, 3. 

Since we have not any information about the behavior of G and A, we "suggest" the following modification. 



B. Modified matter collineations equations 

The problem arises when we are considering a model with G and A variable, so we suggest to change the above 
procedure i.e. the standard one by 

iy(^T,,)=0. (148) 

As in the above case it is easily outlined the resulting matter collineations eqs. so it is a straightforward task to 
integrate them. In this case we have found the following relationships between the quantities: 

^ + ^ = -2^, =^Gp = KoVi-\ (149) 

Gr p Vl 



where Kq is an integration constant, and as above we find that 



V2 = -Vl [Hi-^] X, V3 = -Vl (hi - ^) y, V^ = -Vi (hi-^\ z, (150) 



V = Vidt + Vi\^ ~ Hi] xd^ + Vi(^- H2) ydy + ( ^ - iJg ) zd,, (151) 



obtaining again 
but in this case, it is verified the relationship 

Gp = KoVi-\ (152) 



As in the previous case we may check that its is verified the equation Lho (^^r^T^y^ 



0. 



^{tpG' + tGp' + 2Gp) = 0, ^ 9L + ^ = -l^Gp^t-^, 

G p t 



-px [XX' + tXX" - t [X' f^ =0, ^ tH'i = -Hi X = Xot"' , (154) 

^X^{tpG' + tGp' + 2pG)=0, ^ ^ + - = -- ^Gp«^-^ (155) 



from this equations we also arrive to the obvious conclusions that p = cup with cj G M. 

But unfortunately we have not any information about the behavior of the cosmological constant A for this reason 
we suggest the following modification. 
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C. The complete modification of MC equations. 

In this case we consider that must be satisfied the fohowing equation 

Lv 

that we reinterprets as 

( G{t) , 



(^^T,,-A(t).g,,) =0, (156) 



^T,, =0 = Lv iA{t)g,,) . (157) 



As in the above cases (we are following the same procedure in both cases) we find again that again 

V = Vidt + Vi[^^- xd^ + _ ydy + Vi{^^-H^] (158) 



Therefore checking the relationship 



it is obtained the following results 



Lho (^T,, - A{t)g,, ] - 0. (159) 



SttG 



^(tpG' + tGp' + 2Gp)^-(tA' + 2A), ^ + ^ = _^ ^ « t'^, (160) 

G p t 

p + A] X (XX' + tXX" -t{X'f) =0, tH[ = -Hi ^X = Xot°'\ (161) 



SttG 
SttG 



p + A]y[YY' + tYY" -t{Y'f) =0, ^ tH^ = -H2 ^^ = ^0*"', (162) 



p + A] z (ZZ' + tZZ" - t {Z'f ) = 0, tH'^ = -H3 Z = Zot"^ , (163) 



G' 




G 


P 


tH[ = 


-Hi 






tH',= 




G' 


t = 




P 



^ {tpG' + tGp' + 2pG) = {t A' + 2A) , <=> ^ + = -I ^ Gp ^ , (164) 
c* G p t 

while 

A = Ao^~^ (165) 
that is to say, we have obtained the same results than in the SS section (as it was expected). 

Remark 4 Other possibilities could be explored, for example to calculate the algebra of the matter tensor field T G 
Ti{M) as well as T G Tq(M), note that here we only have studied the case T G r2 (M). Nevertheless and unfortunately 
for these cases we have not been able to find any interesting result i.e. with physical meaning. I would to thank Prof. 
G. Hall for drawing my attention about this fact. 

VII. KINEMATIC SELF-SIMILARITY. 

Kinematic self-similarity has been defined in the context of relativistic fiuid mechanics as an example of incomplete 
similarity (sec for example [54]-[6()]). It should be noted that the introduction of incomplete similarity to general 
relativity is not unique. 

A spacetime is said to be kinematic self-similar if it admits a kinematic self-similar vector V which satisfies the 
conditions 

Cyhij = 25hij, (166) 
CyUi = aui, (167) 

where is the four-velocity of the fiuid and hij — gij + UiUj is the projection tensor, and a and 5 are constants . 
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If 5 ^ 0, the similarity transformation is characterized by the scale-independent ratio a/5, which is referred to 
as the similarity index. If the ratio is unity, V turns out to be a homothetic vector. In the context of kinematic 
self-similarity, homothety is referred to as self-similarity of the first kind. If a = and 5 7^ 0, it is referred to as 
self-similarity of the zeroth kind. If the ratio is not equal to zero or one, it is referred to as self-similarity of the second 
kind. If a 7^ and (5 = 0, it is referred to as self-similarity of the infinite kind, li 5 = a = 0, V turns out to be a 
Killing vector. 

From the Einstein equation (103), we can derive 

CvG,j = SirGCvT.j, (168) 

this equation is called the integrability condition. 

When a perfect fluid is irrotational, i.e., tOij = 0, the Einstein equations and the integrability conditions (168) give 

(a - 5)n,j = 0, (169) 

where TZij is the Ricci tensor on the hypersurface orthogonal to u*. This means that if a solution is kinematic 
self-similar but not homothetic and if the fluid is irrotational, then the hypersurface orthogonal to fluid flow is flat. 

From the physical point of view the detailed study of cosmological models admitting KSS shows that they can 
represent asymptotic states of more general models or, under certain conditions, they are asymptotic to an exact 
homothetic solution [50, 59]. 

Therefore and following the same idea as in the above sections we would like to extend this hypothesis in order to 
find exact solutions to cosmological models with time varying constant. 

Kinematic self-similarity are characterized by the equations (166-167), so in this way it is found that the vector 
field V := KSS is: 

KSS = -{at + P)dt + fixd, + hvdy + hzd,, (170) 

where 

h^(^+{at + P)^, f^=(^+{at + P)^^, /3= |^^ + («t + /3)|^ . (171) 
As in the case of the homothetic vector field in this case it is necessary to satisfy the following ODE 



a—-{at + P) — + {at + P){ — ] =0, (172) 



arriving to the same conclusion i.e 

= 

{at + P) 



Hi = rr^' (173) 



and therefore the solution follows a power law and hence 



X^Xo[t+^) , Y = Yo(t+^) , Z = Zo(t+^] , (174) 



a 

with Xq, Yq, Zq are integrating constants and (ai)i=i ^ We may also check that the KSS vector field must satisfies 
the relationship 

[KSS,^,] = C^jik, (175) 

where ^ is a Killing vector field, in this case = di. This relationship implies that the following equations must be 
satisfy: 



5 + {at + /3) = 0, (176) 



where Hi = ^ respectively, i = 1,2,3. 
In this way we find that 



H ={ai+a2+a2)[t+^) = A (t + 7) \ g = 1 ("1 j _ 1 = _ _ 1, (177) 
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whereA= ^ = 7- 

Therefore, we may calculate the energy density from the conservation equation i.e. 

P = Po(t + 7)-'"'''^^, (178) 
so, following the same procedure as in the above sections, we find that 

A 2 Sirpo 



A' 



(C (t + -iuj + l)AG {t + , (179) 



■2 2A 



with A = Q'iQ2 + 0:30:1 + 02O3, so constant G has the following behavior 

p G 

The behavior of G is the following one 



'0 (t^ + l)A 



as we can see it is verified the relationship Gp « + 7) ^ , as it is expected. 



{growing if {uj + 1)A > 2 
constant if (w + 1)A = 2 . (181) 
decreasing if {lo + 1)A < 2 

with Go > 0. 

Now, we next to calculate the quantity A, therefore from eq. (179) it is found that 



A = Aoit + j)-\ Ao^^{l-j^yfjjjj, (182) 



A 



in this way it is observed that 



Ao > ^ 2 < (w + 1)A 
A = Ao(t + 7)"', { Ao = 0^2 = {lu + 1)A . (183) 

Ao < <^ 2 > (w + 1)A 



The shear behaves as 



(184) 



As in the previous case, the SS solution, we suggest a way to calculate the coefhcients (oi). They must to satisfy 
the field equations i.e. they have to be solution of the following system of equations 

a2 {02 -1) + 03 (as - 1) + a3a2 = A (^—^—^ ' (185) 

ai (ai - 1) + as {as - 1) + a^ai = A (—-^^ ' (186) 

a2 {02 - 1) + ai (ai - 1) + aia2 = A > (187) 

where A = aia2 + o^ai + q;2Q!3, and A ~ ai + a2 + 03. 

So we have the following solutions for this system of equations: 

Qi = a2 = 03, (188) 

Ql = 1 - Q!2 - 03, (189) 



as it is observed we have found the same behavior as in the SS case, but in this case the solution is nonsingular. 
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We would like to stress that a similar solution is already known since 1946 by Narlikar and Karmarkar (see [67]). 
They found the following solution 

ds^ = dt^ - {kt + l)Pdx^ - {kt + Ifdy^ - {kt + Ifdz^, (190) 

where (p, q, r) must satisfy the following relationships: p + q + r = 2, pq + qr + rp ^ 0. 

Before end we would like to make a little comment about the Kasner like solutions. If a solution of (185-187) verifies 
the relationships 

3 3 

^a,2=^a, = l, (191) 

i i 

i.e. they are Kasner's type (see [61], [62] and in particular [44]), then this means that A = aia2 + a^ai + 0203 = 0, 
which brings us to get the following result 



ai = - ( 1 - as - V 1 + 2^3 - 3a2 ) < 0, Vag G (0, 1) , (192) 



a2 = -[l-a3 + Jl + 2a3-3al]> 0, Vag G (0, 1) , (193) 



as = as, (194) 

we think that this class of solutions are unphysical and have a pathological curvature behavior as it is shown bellow. 
We will show that for this class of solutions we have the following curvature behavior 

The non-zero components of the Weyl tensor are: 



C1212 = -«2a3«' {at + ^)-(i+"3+A) ^ ^^^^^ ^ -aiasa^ {at + /3)-(i+"3-hA) ^ ^^^^^ ^ ^^^2 ^ ^)-2(i-"3) 

,+A-l 

(196) 



(03-1)0302 2a3 n OxOsC?, , ^ ^.-(l-aa + A) ^ 0203"^, , _^ (ca+A-l) 

62323 = 5 (ot-t-/^) , 62424 = 5 — (ot-1-/?) ^ 63434= 5 — (ot-l-/?)^ 



where A = \J\ + 2o3 — 3o|. The non-zero components of the electric part of the Weyl tensor are: 

i?22 = {at + /3)-(-+A+i) ^ ^ a^ ^ ^)-(i-..3-A) ^ {l^as) asa- ^ 2(i-.3) 

c c c 

(197) 

and the last invariant has the following behavior 

_ 1- 

c4 {at + P) 

while the gravitational entropy behaves as 



^3 = ''""y.^'-Ti (198) 



(199) 



since I2 = 0. Note that if we take into account another definition for = /3//1, (see [27]-[28]) then we get ^ 00, 
since /i 7^ 0. 

Furthermore, as we can see, if A = 0, then from eqs. (180 and 182) we get 

G = 0, A = 0, (200) 

as it is expected for this class of solutions (vacuum solutions) so they are not interested for us. Nevertheless relaxing 
the condition ^oi^ = \, to our result i.e. ^oi^ <\, we are able to obtain solutions whit (o-;) > 0, Vi, and G ^ 0, A ^ 0. 
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VIII. CONCLUSIONS. 

We have shown how to attack a perfect fluid Bianchi I with G and A variable under the condition divT = 0. In 
the appendix B we will show how to modify these tactics to study this class of models but relaxing the condition 
divT = 0, i.e. without considering the condition divT = and taking into account only the hypothesis of SS. 

With the first of the exposed tactics, i.e. the Lie group one, we have solved the field equations, solving only one 
ODE, eq. (29), studying the possible forms that takes G{t) in order to make eq. (29) integrablc. We have started 
imposing a particular symmetry, X = {at + e)dt + bdp^ to study all the possible symmetries would result a very tedious 
work. 

In this way we have obtained three exact solutions (well actually only two) in function of the behavior of G{t). We 
have seen that the scaling symmetry, X = atdt + bdp, brings us to get the already known flat FRW solution since 
we have arrive to the conclusion that cr = 0, i.e. the shear vanish, and therefore we have rejected this solution since 
we are only interested in solutions that verify the condition a ^ 0. The second of the obtained solutions, exponential 
behavior for G{t), brings us to rule it out, since this solution is only possible if p = 0. The last of the obtained 
solutions, the third one, which is quite similar to the obtained one by Kalligas et al (see [13]), also has been ruled 
it out in spite of seeming with physical meaning, a ^ 0, etc...., since when we calculate the numerical values of the 
exponents of the scale factors (a^)^^]^ , we shown that the only possible solution is the flat FRW one but, which is 
more incredible, with G = const, and A vanishing. This has been a really surprising result, since we think that the 
followed tactic, i.e. solving eq. (29) without imposing any assumption ad hoc, brings us to get consistent results, i.e. 
(T 7^ 0, G{t) 7^ const, and of course, A 7^ 0. 

It is clear that the latter solution is quite similar to the obtained one by Kalligas et al, in fact, we have followed 
their method in order to arrive to the same equation, but all the time, we have tried to avoid to assume ad hoc 
any particular behavior for any of the quantities. Instead of following this way, we have preferred to deduce from a 
symmetry principle the possible behaviors for the function G. Nevertheless Kalligas et at never arrive until the last 
consequences in their calculations, since they did not try to find the possible values for the exponents for the scale 
factors icti)^^i, for this reason their result looks with physical meaning. 

In appendix A we will show by solving the second order differential equation (27), using the same procedure as the 
exposed one in section IV, that it is obtained, at least in order of magnitude, the same behavior as the obtained one 
in section IV. Nevertheless, since we have two constrains, which means that we have two integration constants, and 
therefore in principle they are unknown we are introducing more uncertain in our solutions. Remember that in section 
IV we had only one constrain, Gq, and this fact allows us to arrive to a complete solution for each of the quantities. 
In this way, and following the same procedure, in this case we are not able to rule the solutions out, except in the case 
of exponential behavior, since we get solutions with a ^ 0, although as we already know this is only a mathematical 
drawback. 

Nevertheless and knowing that this method has this kind of drawbacks we have preferred to show both methods in 
other to show that at least in order of magnitude, both methods arrive to the same conclusion i.e. we get the same 
behavior. In a forthcoming paper we study a more complicated model (from the mathematical point of view) which 
is a Bianchi I within the framework of variable speed of light (VSL). In this case will be more useful to study the 
second order ODE instead of the resulting third order ODE, which is really complicate of studying from the point of 
view of the Lie method since this equation will have four unknowns. Therefore, as we already know, we get the same 
order of magnitude studying the second order ODE instead of the third order ODE (other question will be how to 
solve the problem of the integrating constants, but this will be other history). 

At the same time we have shown that it is not necessary to make any ad hoc assumption or to take into account 
any previous hypothesis or considering any hypothetical behavior for any quantity since all these hypotheses could be 
deduced from the symmetry principles, as for example using the Lie group methods or studying the model from the 
point of view of the geometrical symmetries i.e. SS etc... 

With regard to the other tactics employed to study the field equations, i.e. SS, MC and KSS, we have shown that 
both tactics are quite similar and that they bring us to get really similar results, actually as we already know, with 
the SS and the MC we get the same results. 

We have shown that the solution obtained with the SS hypothesis is also quite similar to the obtained one using 
the Lie method under the scale symmetry, except than here, we get the important result cr 7^ 0. This solution also 
is valid for all equation of state i.e. Vw S ( — 1, 1), which enlarge the possibilities for this kind of solutions as we will 
see in appendix B. Furthermore, the exponents must satisfy the following relationships '^cti = 1, and "^olI < 1. 
We would like to point out that if cj = 1 , then we regain the classical solution where G behaves as a true constant 
while A vanish. In appendix B we have tried to show that this class of solutions are consistent with previous results 
obtained by other authors. Nevertheless other authors that have studied this kind of models (under the SS as well as 
KSS hypothesis, but with constants as true constants and therefore with A vanishing) have arrived to the conclusion 
that these solutions must satisfy the relationships X^^^j — 1' ^^'^ '12 '^i — 1- We think that this class of solutions are 
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unphysical since necessarily one of the scale factors must be a decreasing time function (maybe such class of solutions 
would have any interest in the study of singularities). Furthermore such conditions are quite restrictive and impose 
strong restrictions to the curvature tensors, for example, the model is Ricci flat which means that I2 = 0, so the 
gravitational entropy is infinite. 

With regard to the behavior of the "constants" G and A we would like to stress that we have arrive to some 
surprising results since Vcj G ( — 1,1), Gisa decreasing function on time while A is also a decreasing time function but 
negative i.e. A < 0. In the same framework, for the case of the flat FRW, always G is growing. The main diflfcrcnce 
between the SS and the KSS solution is that the KSS one is nonsingular while the SS one is singular. 

We furthermore have pointed out, as it is well known, that if the ST is SS then there is a vector field, V € X{M) 
that satisfies the equation £v9 = 2(7, then such vector field must satisfy the equation £yT = 0, i.e. a homothetic 
vector field is also a MC vector field. Modifying in an appropriate way the MC equations we have been able to find 
the same relationships as in the case of the SS solution. Therefore we have shown that this tactic would be very useful 
in the study of more complicated models as for example the viscous ones. 

To end, we would like to comment the obtained results in appendix B. In this appendix we have found four exact 
solutions for different Bianchi I models under the SS hypothesis. In the first of them, we study the standard Bianchi I 
model i.e. which where G ~ const, and A = 0. We have found again the solution already obtained for many authors, 
but with the restrictions ^Q!i = l,^Q;f<liffw = l. 

In the other studied cases we consider the possibility of one of the constants vary as well as both vary at the same 
time but with the condition divT ^ 0. Therefore, in the model with only G time-varying we find again, as in the 
previous case, that it is only possible if are verified the same conditions i.e. ^ai = l,^af<liffa; = l. Ifcj^^l, 
then the models collapses to the flat FRW one. Nevertheless in the case where A vary or both constants vary we show 
that such possibilities are possible iff ^ = 1 , ^ < 1 and uj ^ 1. We find that G is a decreasing time function 
on time while A is a negative decreasing time function. We have shown how to regain, in a trivial way, the condition 
divT = 0. 
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APPENDIX A: LIE METHOD. 

As we said in the introduction, we would like to compare some different techniques in order to be sure which of 
them is better when one is studying more complicate models, as for example a Bianchi I with c-var. For this purpose, 
in this section, we study the set of solutions for the eq. 

• 2 

p^Ki^ + K2Gp^- - KsAp, (Al) 
P 
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where 



A'i = -— , ^ X3 = 3(l+u.)c2. (A2) 



1 

note that we have the very simple case A'2 = = 1, i.e 



,2 



P^^^-ec^Ap. (A3) 

Observe that we have a second order ODE with three unknowns, which looks simpler than the previous case (see 
section IV), but actually, as we will show, this tactic has many drawbacks, for example the obtained solutions depend 
on many integrating constants and hence it is quite difficult to "predict" their behavior or to get rid of some of them. 

Following the previous procedure we get: 

Ki^P + Pipp = 0, (A4) 
Ki iv - PVp) + iVpp - 26p) = 0, (A5) 
Sp^Cp (KzA - K2pG) - 2R\prjt + 2p\p - ittp" = 0, (A6) 
77p2 {K3A - K2pG) + 2p^^t - K2pG) - pSp (K^A - K^pG) + Cp' {K^A' - K^pG') + p\t = 0, (A7) 

where from (A7) wc get the following constrains if we impose the scaling symmetry X = {at + e) dt + bpdp, (as in the 
above case) 



A' 2a 



A at + e 
G' _ b+2a 
G" ^ ^ at + e' 



(A8) 
(A9) 



Therefore we go next to study the set of the possible solutions. As in section IV, we may follow the standard 
procedure but we would like to point out (to stress) that in this case we have two integration constant Aq and Go 
instead of only one as above so we are introducing more uncertain in this approach. Therefore we find again three 
solutions, the scaling and the full symmetry and the exponential one. The scaling solution has the same behavior 
as the obtained one in section (IV), at least in order of magnitude, but we are not able to rule it out since depends 
of a lot of numerical constants so for example we do not arrive to the conclusion that cr = 0, as above. With the 
regard to the full symmetry we have the same history, we find that this solution has the same order of magnitude and 
that it is nonsingular but depends of many integration as well as numerical constants so the obtained solution is very 
imprecise although it looks with physical meaning. To end and as in the above solution we are also able to get rid of 
the exponential solution, in this case A = Ag = const, and therefore po = 0. We may see, for example, how works all 
this procedure in the case of the scaling symmetry because the rest of the cases are exactly the same as the exposed 
ones in section IV. 



1. Scaling symmetry 



In this case we have X — atdt + bpdp and hence 

A' 2 



^ ^ =^ A = Ao^-^ AoGR, (AlO) 
G at 



G = Gat-'^--, G0GIR+, (All) 



i.e. we are "assuming" that Aq is a real constant (but at this point, we do not know which is its sign) and Go is a 
positive real constant, while 

dt dp b ^1 

-^j- =^ P^Pot^, PoeM+, (A12) 
at bp 

making the assumption, ah < with 6 < 0, and where po is a positive real constant. As we can see, it is verified the 
relationship Gp = i.e. the Mach relationship for the inertia. 
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Therefore 



p xt" is decreasing, G ~ t 



growing if — b > 2a 
constant if — b = 2a 
decreasing if — b < 2a 



but at this point wc have not any information about the sing of Aq. 

Since p must verifies eq. (AlO) with G{t) and A{t) given by eqs. (AlO-All), we find poj 



Po 



c2 {l+u;) + b^^ Sc^Aoa^ (1 + tof^ 



GoUna^ [l + Luy (1 - 
with the only assumption uj G ( — 1,1) ,i.e. lu =/= —1, and cj ^ 1 

Po > <^ b^ > -ab (1 + cj) + Sc^Aoa^ (1 + ujf , 
with ab < 0. Furthermore, p must verify eq. (14) so 



A' 



and hence 



inG' 



5\ 47r 

-Ao = ( 2 + - -^Gopo =^ pa 



oAqC* 



3a (1 + Lj) V ia{l+uj)J \a{l-uj) + 2a + b 



47rGo (2a + b) ' 
2a + & 



(A13) 



(A14) 



(A15) 



(A16) 



(A17) 



Remark 5 As it is observed we are obtaining the same order of magnitude for each quantity, but in this case, we 
have less information about the behavior of the numerical constants, since they depend on more integrating constants. 



As in the above section, with regard to H we find that 

Kt"' ^ Rot-^/<^+'^\ = 

as in the previous cases. 

The shear has the foUowing behavior. 



3 ^ 

Eai = a = — — - — -. 
a[l+uj) ■ 



(A18) 



1 



3c2 [i + ujy \p 



£1 



UJ 



An- 



1 



3c2 (1 + uj) 



t ^a^t , 



(A19) 



therefore we have the foUowing possibihtics Vu; G (—1, 1) . Note that if < — 1 (phantom case) then we need to make 
other considerations. Since > 0, note that t- < 0, then we may suppose that Aq > and then 



6 + 2a < =^ 
6 + 2a = =^ 
6 + 2a > =^ 

but if Aq < then we have the next possibihtics 

62 1 



Ao 

■2 



1 



3c2 {l+oj) 
= 0, for eq.(Al7) 

Ao > 



2 ^1 + 2-^) >0' 



3c2 (1 +w) 



a 



3c^l 



l + 2-j > |Ao| 



(A20) 
(A21) 
(A22) 

(A23) 



with 6 + 2a < 0, otherwise the solutions lack of any mathematical meaning and hence of physical one. In this way we 
note that we have loss of information since we have two integration constant. Go and Ao while in the previous case we 
were able to determine perfectly the behavior of each quantity and in this case we have less information about their 
behavior although both cases are quite similar at least in order of magnitude. 



Remark 6 Therefore, we may say that this tactic is also valid, but has the strong drawback of giving us less infor- 
mation about the obtained solution. 
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APPENDIX B: ON SELF-SIMILAR SOLUTIONS. 

Ill this appendix we are going to study several Bianchi I models and we will show how it is possible to find exact 
solutions to the field equations (without the condition divT = 0) under the hypothesis of SS. 

The time derivatives of G and A are related by the Bianchi identities i.e. eq. (5) that in this case collapses to the 
following one: 

P + p(l+^)(^f + 7 + 1^ =/W, (Bl) 

where f{t) is a function that depends on time and controls the time variation of the constant G or/and A. If G = const. 
and A vanish then fit) = 0, so the model collapses to the standard one. This idea was pointed out by Rastal (see 
[68]) and improved (in the theoretical framework of time varying constants) by Harko and Mak (see [09]). 
Therefore the resulting field equations are (9-12) together to the new one 

p + p{l+u:)H = f{t), (B2) 

with 

= (^1 + I + 1^ = 3| and i?-^ = XYZ, (B3) 

i.e. we are following the same notation of the previous sections. 
We remember that the homothetic vector field is: 



V = tdt + \ l-t^\xd^ + (l-*^) ydy+ ) (B4) 



so this means that we have obtained the following behavior for the scale factors: 

X = Aoi"! , Y = Fot"" , Z = Zaf^' , (B5) 
with Xq, Yqj are integrating constants and (cti)^^-^ G R. In this way we find that 



1. "Constants" constant. 

In this case we consider f{t) ~ 0, so this means that G = const, and A vanish and therefore we get that from cq. 
(B2) that 

p + p{l+uj)H ^0, ^ p = poi~^"+'^". (B7) 
From the field equations (9) we get that 

Ac^ 2 
SttG [uj + 1) 

where A = 0L\a2 + ol2,oi\ + ol20lj,. 

The shear has the following behavior, ^ 0, as it is observed tr ^ as (a^ — > oij) . As in the previous sections, we 
may calculate the coefficients (a.i) by solving the following system of equations: 

a2 {oL2 - 1) + as (q!3 - 1) + a.za2 = -Au, (B9) 

a.\ (a\ — 1) + as (as — 1) + asai = —Au, (BIO) 

a2 (a2 — 1) + ai (ai — 1) + aia2 = —Au, (BH) 

a(w+l) = 2, (B12) 
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where A = a\a2 + asai + a2a3, and a = cx\ + cxi \- a^- 

So we have the following solutions for this system of equations: 

ai = 1 — Q!2 — as, w = 1, (B13) 
2 

ai = a2 = as = 177 — — rr, (B14) 
3 (tj + 1) 

as it is observed only solution (B13) is interesting for us. The second solution is the usual FRW one, so it is not 
interesting for us (see Einstein&de Sitter ([(iH]) for a; = 0, and Harrison ([f)4]) Vo;). Nevertheless we have found that 
solution (B13) verifies the conditions 

a = ^ a, = 1, ^a\<\, (B15) 

but iff w = 1. (see [43]) while other authors claim that must be satisfies the condition X)*^? = I7 (see [Gi], [62] and 
[65]) and in particular, in this context (see [44]). 

Therefore we have obtained the following behavior for the main quantities: 

i?=p =^ g = 0, (BI6) 

so it is quite difficult to reconcile this model with the observational data. With regard to the energy density we find 
that 

and with regard to the constants (Q:i)i=i '^e have that only obtain a BI solution iff ol\ ~ l — a2 — a3, (where furthermore 
we suppose that a2 ^ as) and that this result only is possible if the equation of state is w = 1, i.e. ultra-stiff matter 
(see [43]). For a review of Bianchi I solutions see for example ([6()]). 



G— variable. 



In this case we are going to consider that only vary "constant" G. This only possible if we take into account the 

G ' 



condition divT ^ and therefore f{t) ~ ^ttP, so eq. (B2) collapses to the following one 



- + ^ = -(1 + ^^)7, =^ pG = t-(i+-)", (B18) 
p G t 

From the field equations (9) we get that 

r'^ A ^ 2 
Gp^ 7 7^, a = - -. (B19) 

The shear has the following behavior, ^ 0, as it is observed cr ^ as (a; — > aj) . As in the previous sections, to 
calculate the coefficients (a^) we need to solve the following system of equations: 

q;2 (a2 - 1) + as (as — 1)4- asa2 = —Auj, (B20) 

ai (ai - 1) + as (as — 1) 4- asai = -Auj, (B21) 

a2 (a2 — 1) 4- ai (ai — 1)4- aia2 — —Au, (B22) 

a(w+l) = 2, (B23) 

where A = aia2 4- asai 4- a2as, and a = ai + a2 + as. So we have the same solution as in the above case. 
Therefore we have obtained the following behavior for the main quantities: 

H=j, =^ (B24) 

and 

G,= ^t- -^=3^(1-3A)1 (B25) 

note that this result is quite similar to the obtained one in the last solution i.e. the obtained one in eq. (B17), but 
we are not able to get a separate behavior for the quantities G and p. 
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3. A variable. 



In this case we consider only the variation of the cosmological constant A, so eq. (B2) yields 

, f X Y Z\ Ac^ X 
P + Pd + c.) - + - + -U-_, (B26) 



and therefore from the field equations (9) we get that 

c2 c' 



, A 2 SnG I 
^ = --7^ - — P' (B27) 



and hence 



Now, we next to calculate the quantity A, from eq. (B27) we get 



in this way it is observed that 



r Ao > ^=> {uj + \)a > 2 
A==Ao^"^ <^ Ao = ^ (cj + 1)q' = 2 . (B30) 
[ Ao < {uj + l)a<2 

The shear has the following behavior, ^ 0, by hypothesis. As in the previous sections, we calculate the coefficients 
(oii) from following system of equations: 



CC2 {o-2 - 1) + as ("3 - 1) + Oiza.2 = A \^ ^ j ' (B31) 
ai (ai 1) + Q!3 (q!3 - 1) + as"! = A — , (B32) 
0.2 {ct2 " 1) + ai (ai - 1) + a\ai = A \ ) , (B33) 



where A = OL\ai + aaai + a2a3, and a = ai + 02 + 03. 

So we have the following solutions for this system of equations: 

ai = a2 = Q!3, (B34) 
ai = 1 - q;2 - Q!3, (B35) 

as it is observed solution (B34) is not interesting for us, since it is unphysical (in this context). Only the second 
solution has physical meaning and it is valid Vw G (—1, 1]. Nevertheless we have found that this solution only verifies 
the first of the condition of the Kasner like solutions i.e. 

a = ai = 1, and of < 1, (B36) 

and it is valid e ( — 1,1]. 

Therefore we have obtained the following behavior for the main quantities: 

H=j, =^ g = 0, (B37) 



while with regard to the energy density we find that 
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so, if a; < — 1 =^ p is negative (phantom cosmologies), for the rest of the values of i.e. uj G (—1, 1], p is a decreasing 
fimction on time. 

The cosmological "constant" behaves as follows 

so we have found that A is a ''''negative decreasing function^^ on time. As we can see this solution is quite different of 
the previous ones, since here we have obtained a solution type Bianchi I V w G (—1, 1) while in the previous ones this 
only happens if a; = 1. Here if = 1 then we regain the first of the studied cases i.e. which one where A vanish and 
G behaves as a true constant. 



4. G&A-variable. 



In this case we are going to consider that both "constants" G and A vary, as in section V, but in this occasion it is 
not verified the additional condition divT = 0, therefore eq. (B2) yields 



SttG 



P + + f + 1 + 1 



-A - —Gp, 



so from the field equations (9) and (14) we get that 



and taking into account the eq. 



we get 



p + p{\+uj)H = - 



Ac^ G 



SttG G' 



2A 



Stt 

72" 



-Gp{l + u)- 



and hence 



Gp 



A 



-r 



47r (cj + l)a 



(B40) 

(B41) 

(B42) 
(B43) 
(B44) 



as we can see it is verified the relationship Gp « t^^ , as it is expected. In fact it is impossible to separate both functions 
(to get the behavior of both functions independently), to do that we need to impose a condition, but precisely we are 
trying to avoid such way. 



Now taking into account again eq. (9) we get 



Ao = 4 (1 



{lo + l)a) ' 



in this way it is observed that 



A = Aor 



Ao > <F=^ [lo + 1)q; > 1 
Ao = {oj + l)a = 1 
Ao < [lo + 1)q' < 1 



(B45) 



(B46) 



The shear behaves (see eq. (20)) as follows: ct^ ^ 0, by hypothesis. In order to find the value of constants (a^) , we 
make that them verify the field eqs. so in this case we get the following system of eqs.: 



Q!2 ("2 - 1) + "3 (q!3 - 1) + a3a2 = A 
Oil (Q!i " 1) + 0^3 (a3 — 1) + a^ai ~ A 
a2 (ot2 — 1) + ai (a\ — 1) + q;iq;2 = A 



a 




2 




a 




a 




2 




a 




a 




2 



(B47) 
(B48) 
(B49) 
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where A = aia2 + ctsai + a2a3, and a = ai + a2 + (X3. Therefore we obtain the same solution as in the last studied 
case finding in this way the following behavior for the main quantities: 



and with regard to the product Gp we get 



H=^: =^ q = 0, (B50) 



but we cannot say anything more. The cosmological "constant" behaves as follows 

^=^ot-^ Ao-{i:;S,""v:e(-M) ' 

so we have found that A is a negative decreasing function on time. 

In order to try to find a separate behavior for the functions p and G, we may suppose that 

P = Poi"", G = Goi"-2^ =^ Gp^^-^t-^ ^Kt-\ (B53) 

47r yijj + 1) 

with a e M^, i.e. for example we may choice 

47r/9o (t^ + 1 ) Po 

therefore, it is verified the field eq. (B4G) for all the possible values of a, but if a = (w + 1) then we regain the 
condition divT = as well as /(t) ~ 0, i.e. 

dwT^p + p{l+uj)H^Q^-^-^p^f{t). (B55) 



